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Abstract. The Schwinger equations of QED are rewritten in three different ways as integral equations 
involving functional derivatives, which are called weak field, strong field, and SCF quantum electrodynam- 
ics. The perturbative solutions of these equations are given in terms of appropriate Feynman diagrams. 
The Green function that is used as an electron propagator in each case is discussed in detail. The general 
renormalization rules for each of the three equations are provided both in a non perturbative way (Dyson 
relations) and for Feynman diagrams. 

PACS. 12.20.-m Quantum electrodynamics - ll.10.Gh Renormalization 
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This paper is a step towards the calculation of photon 
and electron spectroscopies of matter based on quantum 
electrodynamics. 

Starting from quantum electrodynamics, which is a 
most accurate and successful theory, garanties that the ba- 
sis of the calculation is sound. Although such an approach 
may look too true to be beautiful, it seems adequate be- 
cause of experimental and theoretical reasons. From the 
experimental point of view, spectroscopy has made huge 
advances and is now able to measure tiny effects. From the 
theoretical point of view, relativistic methods have proved 
very powerful, even for problems that do no look obvi- 
ously relativistic. For instance, relativistic density func- 
tional theory Jlj] , ^ does not meet the nonuniqueness prob- 
lem of spin-density functional theory ||. Moreover, rela- 
tivistic quantum field theory is a safe framework to go 
beyond the present LDA methods, especially when deal- 
ing with excitations. 

There are many presentations of QED. We first have to 
choose the version which is most convenient for solid-state 
physics. The main three formalisms used in textbooks are 
the quantum field, path integrals and Schrodinger rep- 
resentations 0. The quantum field and Schrodinger ap- 
proaches manipulate operator-valued distributions g and 
the path integral approach uses ill-defined measures. There- 
fore, we prefer to use Schwinger's idea Ref.Q of working 
directly with Green functions which are given as solutions 
of equations involving functional derivatives. 

This approach has the mathematical advantage of ma- 
nipulating standard distributions, and the physical advan- 
tage of using only measurable quantities, or the more mea- 
surable ones. 



The standard approach to QED comes clearly from 
particle physics, where the S matrix is most useful, and 
where in and out states of the scattering experiments 
are well defined. In solid state physics, measurements are 
usually based on a different principle. The spectroscopist 
shines on the sample a beam of electrons or photons com- 
ing from a classical source. By classical we mean that the 
source is not influenced by the system being measured. 
After its interaction with the sample, another beam of 
electrons or photons is measured. Many spectroscopies 
can be described within this framework: photoemission, 
electron scattering, photon scattering, inelastic scattering, 
BIS, LEED, RHEED, x-ray absorption, UV/visible spec- 
troscopies, etc. We try to stick as much as possible to this 
experimental point of view. 

We describe the photon field with the electromagnetic 
potential (which is not directly measurable but can be 
calculated in a given gauge from the measurable elec- 
tric and magnetic fields) and the photon Green function 
(from which the photon energy density can be computed) . 
We describe the electron field with the one-particle Green 
function. This is also not directly measurable, but many 
diagonal matrix elements of the Green function can be 
measured (e.g. the electronic charge and current densi- 
ties). 

For applications to the spectroscopy of matter, we need 
to formulate QED with an external field. This external 
field is made of the potential due to the nuclei, to mag- 
netic or electric fields applied to the sample and to ex- 
ternal light sources. After this introduction, we start with 
the definition of the notation used for the Green functions 
and the QED Lagrangian. Then, we derive the Schwinger 
equations for the electromagnetic potential and the elec- 
tron propagator. These equations are solved iteratively us- 
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ing three different methods (weak field, strong field, self- 
consistent field). The boundary conditions, such as the 
number of electrons in the system, are determined by the 
unperturbed Green function which is discussed in detail. 
The proper definition of the current is established through 
the use of a bilocal operator. Then we describe the renor- 
malization of QED in an external field. The non perturba- 
tive Dyson relations are given to express the relation be- 
tween bare and renormalized Green functions. The renor- 
malization rules for Feynman diagrams are discussed. 

The following assumptions are made in the present pa- 
per: the external field is weak enough not to create charges 
(this is true for all stable atomic nuclei 0) and the ex- 
ternal field is zero at infinity (we exclude constant electric 
and magnetic field, for which specialized monographs are 
available ||,[[§). Moreover, we do not consider IR diver- 
gences. 



2 Notation 

In this section, we specify the notation that is used in 
the paper. The charge of the electron is e = — \e\. The 
pseudo-metric tensor is 



where 



9 = 



(\ 

0-10 
0-10 
\0 -1, 



We choose the standard gamma matrices (pQ] p. 693) 
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The Pauli matrices ex-? are defined by 
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The charge conjugation and mass reversal matrices are, 
respectively 
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-ia 2 \ , (0 1 

- lf ? 2 ' and 75= U0 



They satisfy the following identities 
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2.1 Free propagators 

The free photon Green function D^(x, y) is defined by its 
Fourier transform 



Dl v {x,y) 



<W ex V [-iq.(x-y)]D°( q ), 



KM 



1 



—{(9^ - (i - i/£)<?m<W9 2 )- 



e (g 2 +ie)' 

It is a solution of the equation 

e (D<r - (1 - Od»d v )D vX {x,y) = S\S(x - y). 

In the last equation the D'Alembertian □ and the deriva- 
tives d^d u act on variable x. 

The free electron propagator is 



S°(x,y) = 



dq 



where 

S°(q) = (hcj ■ q - mc 2 + it)" 1 . 

It is a solution of the equation 

(ifkr/ ■ d — mc 2 )S°(x, y) — 8{x — y). 

Analytic expressions for these (and other) propagators 
are given in |nj and [jl2| section 2.3. Notice that these 
Green functions arc neither advanced nor retarded. Their 
physical meaning is discussed in jl3j, (lj] and (l5) p. 57. 

2.2 The QED Lagrangian in SI units 

To compare our results with those of nonrelativistic many- 
body theory, it will be useful to write the QED Lagrangian 
in SI units p|. 

We have xr = ct, we define the 4-current J a — (p, J/c) , 
and the 4-potential A a = (V,cA). 

The field-strength tensor is (j|10| p.8) 

-E x -E v —E z 

F a/3 _ qq a 13 _ a/3 Aa _ I -tf3 z cB y 

1 /•;'' cB z ell' 
E z -cBy cB x 

The Maxwell equations are (jl0| p. 9) 
8 aF »P = lj(\ 

The photon Lagrangian is (JIo| p. 12) 



£ 7 = -|f q/3 F^ = |(|E| 2 -c 2 |B| 2 ). 
The electron Lagrangian is 



(2tt 



C e = ihctpj ■ dip — mc 2 'tp'ip, 
where 7 • d = 7^9^. The interaction Lagrangian is 
Ci = — J • A = —etpj ■ Atp. 
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The current is 



JJx) = eip(x)^/a^(x). 



The total Lagrangian in SI units, including the gauge term 
is 

C = - *jF aP F af) + tp(ihey ■ d— ej- A— mc 2 )^ 



2.3 Dimensions 

The following table gives the SI units of the quantities used 
in the paper, in a space-time with dimensions d (d = 4). 



3 Derivation of the functional equations 

The Schwinger equations were presented in Ref.jq], and 
various derivations of them are available p^| , Jul (p. 416- 
32) and (p. 475-81). Our derivation follows Ref.|o|. 

We define a generating function Z = Z(j,r),fj), where 
the photon source j^ix) and the anticommuting electron 
sources rj(x) and fj(x) enter the total action as |L0| 



(7 = 1 



dx j fl (x)A tJ '(x) + / dxi](x)ip(x) + ip(x)r](x), 



where / = J dxC(x). The minus sign before j^,(x) in the 
definition of a was chosen so that j^(x) is a standard elec- 
tromagnetic current. It is the opposite of the convention 
used in Ref.@. 

The generating function can be written, up to a nor- 
malization factor, as a path integral (pfl p. 476) 



Z(j,V,v)= / £>(A^)e iCT/ftc 



(1) 



or as the mean value of an operator (JhJ p. 210, 261) 

z {j,fj,T]) = (Texp (i(a -/)/&)). 

It is also possible to write Z(j,fj,r]) in terms of a 
ZotiiViV) without interaction Jl(J p. 445 and |l§ p. 246. 
Anyway, we do not calculate Z, we only need it to calcu- 
late mean values of Heisenberg operators. For instance 

Z(A^(x)) = [ V(A^,^)A^x)e^ hc . 



From the definition of a, A^x) can be written as a func- 
tional derivative with respect to j^{x): 



Z(A fl (x))=- / V{A^,$)- 



8a 



8j»(x) 



7 /he 



= ihc- 



SZ 



8j»(x) 
Similarly, we shall use 

Z{TMx)Mv)) = 



(2) 



5 2 Z 



Sr)P(y)5fj a (x) 



and 



ZiTA^A^y)) = -(he) 2 



S 2 Z 



8j»(y)8j x (x)' 

Finally, the following property will be essential 

M*)> = $(x)) = 0. 

These equations are derived from the fact that the ground 
state of the system is an eigenstate of the charge. If Q is 
the charge operator, then [Q,ip(x)] = —tp(x) (jlQ] p. 147), 
thus 

(V(a:)) = (G\i>{x)\G) = -(G\Qyj(x)\G) + (G\^(x)Q\G). 

The ground state is an eigenstate of Q with eigenvalue N, 
thus 

(G\i/){x)\G) = -N{G\i)(x)\G) + N(G\if)(x)\G) = 0. 

More physically, no anticommuting operator is measur- 
able, only products of an even number of anticommuting 
operators can be observed. It may be stressed that the 
mean values of ip(x) and ip(x) are zero only when the (un- 
physical) electron sources r](x) and fj(x) are set to zero. 



3.1 The photon equation 

To derive the Schwinger equations we use the fact that 
the derivative of an integral is zero (assuming that the 
integrand is zero at infinity): 



V(A,i>,iP)^^ = 0. 



Thus 



81 



,8A^(x) 
A direct calculation leads to 
81 



8A»(x) 



^ia /he 



(3) 



8A»(x) 



e {Ug^-{\-Z)d vl Bu)A v (x) 
-e^(x)7 A ,V'(x). 



(4) 
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According to the method of generating functions, a 
factor A u (x) in the integral can be replaced by a func- 
tional derivative of Z with respect to j v (x) (up to a factor 
ihc). This is what we did for Eq.(||). In the case of an- 
ticommuting variables such as ip(x) we must be a little 
more careful. Functional derivative with respect to an an- 
ticommuting source is very similar to that with respect 
to a function. The difference can be summarized in the 
identity 

5 (AB)-^ A B + ( l) lAl A SB 
Sr] Sr) Sr) ' 

where A is a product of \A\ anticommuting variables. For 
instance, \A^\ = 0, = 1, \f}\ =1, \tp^ip\ = 2, etc. 

Thus, each factor ip(x) (resp. ij}{x)) in the path integral 
is replaced by a functional derivative of Z with respect to 
fj{x) (resp. f)(x)). Paying attention to the signs and the 
factors hc/i we can rewrite Eq.(||) as 



SI 



hcS hcS —hcS 



8A^(x) V iSj ' iSrf iSrj ) ^ 
Therefore, Eqs.(|) and (|) in Eq.(|) yield 

— hc5Z 



(5) 



eo(Og^ - (1 - Qd^dv) 



iSj v (x) 
hcSZ 
i8rj s (x) iSfj s ' (x) 



—hcS 
-e — - , 7, 



Zjfj,{x). 



The value of the vector potential and the electron wave- 
functions are obtained by |l(J 



he 1 SZ 

1 Z 5j fJi (x)' 

he 1 5Z 

i Z Sf) s {x) 1 

T . . he 1 SZ 
tp s {x) = ~— . 

2 Z QT] s (x) 



(6) 
(7) 
(8) 



To simplify the notation, we have written A^x), ip s (x) 
and 4> s (x) for (A M (a;)), (ip s (x)} and (ip~ s (x)). 
Thus, we obtain 

Ze (Og^-(l-Z)d^)A»(x) 

-eiTic j s / Zip s > (x) = Zj^(x). 

We define now the electron Green function by 

Sss ,(x,y)=ihe 6H ° gZ 

5r] s ' (y)Si] s (x) 

Sip s (x) _ Sip s > (y) 
5r) s ' (y) 5fj s (x) 

Thus, the equation becomes 

eo(ng tlv -{l-Z)d lt d v )A v (x) 
-e-tp{x)^^ij}{x) + iehctr[y fJ ,S(x, x)\ = j^{x). 



(9) 



When the external electron sources rj and fj are set to 
zero, we showed that ip(x) = ip(x) = 0, and we obtain our 
first basic equation 

e {u 9llv -(l-^d v )A v {x) 

+ieHctr[j fl S(x,x)] = j^x). (10) 

It is a bit clumsy to use an electron source rj(x) just 
to conclude that no such source exists which leads to the 
cancelation of ip(x) and ip( x )- The way out of this difficulty 
is to use Rochev's bilocal source T)(x, y) J2(J. 

cr = 7- / dxj^A^ix) + / dxdyi)(x)r](x,y)i)(y), 



where r)(x,y) is now a physically reasonable source of 
electron-positron pairs. Such a source would lead imme- 
diately to Eq.([To|). We used the more standard electron 
sources to follow the textbook derivations jlJJ. 

Equation ( |l0| ) means physically that an induced cur- 
rent — iehctr[ r y^S{x 1 x)] must be added to the external cur- 
rent jn(x) as a source of electromagnetic potential. 



3.2 The electron equation 

The second equation is obtained by varying ip(x). Follow- 
ing pj p. 478 we obtain 



51 r -hcS hcS -hcS \ = 
Sip{x) V iSj ' iSf) ' i5r/ J 

The functional derivative of the action yields 

51 

— -=-r — = (ihcyy ■ d — mc 2 )ijj{x) — ej ■ A(x)ip(x). (12) 
Sip(x) 

Introducing this into Eq. ( Jl~T| ) , we obtain 

2 hcSZ —HcS hcSZ 

—Zr)(x) — [incy ■ o — me )- 

From Eq.(^) we can write 
he SZ 



iSrj(x) ' ^ iS j fj,(x) iSrj(x) 



i Sr) s (x) 



Ztp s (x). 



Hence 



-r](x) = (iher/ ■ d — mc 2 )il){x) — iehcy^ 
—ej ■ A(x)ip(x). 



8ip(x) 



Finally, we take the functional derivative of this equa- 
tion with respect to T)(y), we use Eq.(^), the fact that 
ip(x) — and we obtain our second basic equation 



S(x — y) — (ihey ■ d — me 2 — ■ A(x))S(x, y) 
SS(x,y) 



-iehcj^- 



(13) 
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Again, the same result can be obtained without elec- 4 Three integral equations 

tron sources 77, fj by using the electron-positron source 

v( x iV) [§0)' j n this section, we derive various integral equations which 

The equation for the electromagnetic vector potential correS pond to the differential equations of the previous 

IS 

section. 

A„(x) = o p (a) - iehc [ dsD° (x, a)ix[fS(8, s)}, (14) 



, /•%•,! , 1 , ,• 1 iii ii 4.1 Weak external potential 

where a^(x) is the external potential created by the ex- r 

ternal current j v (x): 

When the external potential is weak, we can multiply 

a (x) — / dyD°„(x, y)j u (y)- Eq.(^5|) by the free electron Green function S°(x, y). This 

J gives us the following coupled equations 

In terms of the external potential (x) , the functional „ 

equation (@ becomes S(x, y) = S°(x, y) + e dz5°(x, z) 7 • A(z)S(z, y) 

(ihc-/-d-mc 2 -e-/-A(x))S(x,y)=S(x-y) r SS(z,y) 



XQ( r „) + ieh c / dzdz'S {x,z)^ ' D (z',z), 



im Di ^ x) - (i5) 



On the other hand, we can take the functional deriva- 
tive of the action with respect to ift( x ) 



A^x) = a^x) - iehc / dyD° (x, y)tr[y v S(y, y)]. 



The weak field approach to atomic physics was re- 



= ihcd^x)^ + $(x)fa ■ A"{x) + mc 2 ). (16) viewed recently by Eides and ColL @ 

oip[x) 



We denote the external current j /1 (x) by a star * , so 
the external potential a^^x) is denoted by the Feynman 



We can now repeat the calculation that was done start- 
ing from Eq.(h2). This gives us another equation for the 
electron Green function diagram . 

S(x — y) = -ihcd£S(y, x)^^ - S(y, x)(ej ■ A(x) + mc 2 ) 
h SS(y,x) 

le Sj (x) 4.1.1 Feynman diagrams 

The coupled equations generate the following series for the 
3.3 Photon Green function electron Green function. The propagator is oriented from 

. ..... right to left and the electron loops are oriented anticlock- 

For spectroscopic applications, it is useful to know the w j se 

photon Green function 



6j"(y) ' e° -> • ■ = 5° 



In physical terms, the photon Green function gives the lin- 
ear response of the electromagnetic potential to the vari- e 1 
ation 5j v (y) of the external source by 



5A li (x) = I dyD flv (x,y)6j v (y). e 



2 



A 



rr 



To obtain an equation for D^ v {x, y), we solve Eq.(|lO[) 
for A^x). ^3 _^ 



A ll {x)= I dzDl x {x,z)(j x {z)~iehctr[ 1 x S{z,z)] 



TTT 

A functional derivative with respect to j v (y) gives us the + ? + 

equation for D^(x, y) • 1 £ — t ■ . C — 4 

D IMV {x,y)=Dl v {x,y)-iehc f dzD^x, z)tr[ 7 ^ S{ ' ] 



(18) 
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4.2 Strong external potential 



When the external potential is strong, we multiply Eq. ( |15| 
by the Green function in the presence of a^(x): 

S a = (ihcj ■ d — mc 2 — ■ a(xj) 
Now the coupled equations become 



S(x, y) = S a (x, y) + e dzS a (x, z) 7 • A'(z)S(z, y) 



+iehc J dzdz>S a (x,zh/-^}±D%(z',z), 
A'Jx) = -iehc [ dsDl u (x,y)tv[YS(y,y)}, 



In this case, the total potential is ^(a;) = a^(x) + A'(x). 

Strong field QED was reviewed recently in Ref. p2fl . 
Strong field QED can easily accomodate bound states, 
but the nuclear potential is not screened. 



4.2.1 Feynman diagrams 



If s n (resp. a n ) denotes the number of Feynman dia- 
grams for S(x,y) (resp A^(x)) at order e™, the coupled 
equations yield 



•so 

On 



1, a = 0, 



Sn— 1 , 
71-1 



S„ = Sidn-i-1 + (n - l)Sn-2- 



Thus, s 2n = 1,2,10,74,706,8162,110410,... and s 2n +i = 
0. If s(x) is the generating function for the sequence s n , 
then 



For the electron Green functions, the propagator is ori- 
ented from right to left and the electron loops are ori- 
ented anticlockwise, except for the tadpoles, where a loop where 
diagram is half the sum of a clockwise loop and an anti- 
clockwise loop. 



s(x) = ^ s nX n = J 



71=0 



f(x) = £(2n+l)!!a: 2n . 



For similar results, see 



p.467. 



=0= 




+ 




4.3 SCF external potential 

Finally, the most accurate results are obtained when the 
external field is taken as the complete vector potential 
A n (x), which is determined by a self-consistent field pro- 
cedure. Therefore, the potential that will be used in the 
initial Green function is A^(x) instead of a^{x). This cal- 
culation can be found in Refs.[^8| and |l0| p. 479. 

To eliminate the external current j u (y) from Eq.(13) 
we need the functional relation 



SS(x, z) 

WW) 



ds 



SS(x,z) SA^s) 
SA^s) 6j»(y) 



SS(x,z) 

As Jam d ^ v) - 



With this relation, the functional equation for the pho- 
ton Green function Eq.(|l8|) becomes 



(x, y) = D° (x, y) - iehc \ dzdsD° x (x, z) 



For the potential A^ (x) we have the following Feynman 
diagrams 



. ySSfz, z), . . 
oAx'(s) 



(19) 
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and the functional equation for the electron Green func- 
tion Eq.([l^) becomes 

(ifkrf ■ d — mc 2 — ey ■ A(x))S(x, y) — 8{x — y) 

+iehcrf [ d /S(x g 
J oA x {s) 

Then, we multiply Eq.(po|) by the Green function in 
the presence of A fJi (x): 

S A = (ihcy ■ d — mc 2 — e7 • A(x)) , 

and the equation for the electron Green function becomes 

S(x,y) = S A (x,y) 

+iehc J dzdz'S A (x,z)- / J^0^D^(z' 1 z).(21) 

This equation, together with Eq.(|l9|) and Eq.(p^) are a 
complete system of equations for the determination of the 
bare Green functions and the potential. 

As compared to the case of a strong external poten- 
tial, the SCF potential has the advantage of taking into 
account the electrons in the system. Thus, the nuclear po- 
tentials are screened by the electrons in a self-consistent 
way. All the tadpoles of strong field QED have been re- 
summed. 

This formulation of QED is closer to the standard 
methods of solid state physics or quantum chemistry. It 
was used in nuclear physics under the name Hartree QED 
when the current is calculated by a single electron loop 
|| , see also @. 



4.3.1 Electron Green function for self-consistent field 

For the electron Green functions, all electron loops are 
oriented anticlockwise and the propagator is oriented from 
right to left. 




The number of Feynman diagrams at each order was 
calculated in Ref. |25|] (see also pq ). 



For the potential, we have the following diagrams 



e 3 - ~h£> 

Each strong field diagram is the sum of an infinite 
number of weak field diagrams, and each SCF diagram is 
the sum of an infinite number of strong field diagrams. 



5 The initial Green function 

Once the functional equation is know, it remains to specify 
the boundary conditions. For example, we can specify the 
Green function without interaction. In the case of a static 
external field a /i (r) or a static self-consistent field A M (r), 
this Green function is well-known and can be written 

G Q N (x,y) = i9(y a -x () ) £ $ n {*)Mv) 

E n <E F 

- t 8(x°-y°) M*)$n(v), ( 22 ) 

E n >E F 

where Ep is the Fermi energy, and where ip n {x) is the 
solution of the Dirac equation at energy E n . Note that 
G%(x,y) = hcS%(x,y) or G° N (x,y) = hcS$(x,y) where 
S^(x,y) and S^(x,y) are the Feynman electron Green 
functions in the external potential a^{x) or A^{x) (with 
N occupied bound states). In the present section, we want 
to discuss and justify this expression. 

The notation follows standard textbooks on the Dirac 
equation The index N 

means that G° N (x,y) is 
the Green function for a total charge of Ne. The Green 
function Gq(x, y) represents the vacuum in the presence 
of tfy(x) (or A M (x)). This vacuum Green function is given 
by Eq.(p2[) for Ep — 0. The difference between G° N (x,y) 
and Gq(x, y) is 

G%(x,y)-G° (x,y) = i6(y -x ) £ M^)Mv) 

0<E n <E F 

+i9(x°~y°) ^n(x)My), 

0<E n <Ep 

= * ^2 V>n(a#n(y)- 
0<E n <E F 

The Fermi energy Ep is determined by the condition that 
the sum is over TV states: 

E ! = 

0<E„<E F 
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5.1 Early works 

To understand the origin of the Green function given by 
Eq.(p2|), it is useful to describe how it appeared histori- 
cally. 

In the vacuum, Dirac p8j,]2j| assumed that all the 
negative-energy states are filled. This gives the density 
matrix 

ot. 



P (x,y)= E <( x )V>f(y), 



£„<0 

where V'nW is the free solution of the Dirac equation 
for energy E n . However, such a density matrix would be 
highly impractical, since it would lead to an infinite charge 
density. In a strike of genius, Dirac made the following 
observation Q,|2j|. Since 

^V°(x)^ f (y)=<5(x-y), 

the "renormalized" density matrix 
p (x, y) = p(x, y) - ~<5(x - y) 

E„<0 E n >Q 

is much better behaved^ In particular, the charge density 
is 

p (x) = tr[p (x,x)] 

= JE \^)\ 2 -\ E iv^(-)i 2 = o, 



where the last equality is due to charge conjugation sym- 
metry. 

In the presence of a static potential a M (x) or ^(x), 
and without bound states, the analogous density matrix 
is 

p{*,v) = \ E M*)€\y) -\ E MWnfrl 



E 7l <0 



E n >0 



where ip n {x) is now the solution of the Dirac equation with 
a potential a^(x) or ^(x) and energy E n . Of course, the 
charge conjugation symmetry is broken by the external 
potential, but we still have that 

p(x) = tr(p(x,x)) = tr(p(x,x)) - p (x) 

= \ E(l<M*)l 2 -h/tol 2 ) 

E n <0 



\ Ed^WI 2 "!^^)! 2 ) 



E n >0 



1 In fact, Dirac defined po( x , y) as the double of the right- 
hand side. This was immediately corrected by Heisenberg |30|] . 



is less singular. The physical idea developed by Dirac is to 
measure the charge density with respect to charge den- 
sity of the vacuum po(x). The induced charge density 
p(x) ^= in the presence of an external potential was in- 
vestigated in the thirties , @ , |f| , Q , (32| . The po- 



tential corresponding to this induced charge is called the 
Uehling potential Q|. Physically, p(x) ^ is a reaction of 
the vacuum which is polarized by the external potential. 
The induced charge has observable consequences. For in- 
stance, its effect on the energy levels of hydrogen atoms is 
well documented ]22[ ] . Computer programs are available to 
evaluate the Uehling potential for general nuclear charge 
models [Q. Expansions in (Za) n where investigated in 



5.2 Checking the Green function 

In this section, we check that Eq.(|22|) gives the charge 
density proposed by Dirac. In Eq.(|22|) we write 

ipn(x) = exp(-iE n x° /ch)4> n (x-), 
and we define the frequency dependent Green function 



where, on the right hand side, x° = y° + ct. Since the 
external potential is time independent, the Green function 
depends only on t. To determine G^r(x, y; lo) we use the 
distribution identity (Ref.JTf|, p.92) 



9(x° - y°) = i 
and we obtain 



G^(x, y;w ) 



dw exp [— ioj(x° — y )] 



2tt 



lo + it 



E 

E n <E F 

+ E 



^n(x)^„(y) 
lo — E n /H — it 



E n >E F 



lo — E n /H + it 



This is the standard expression for the non-interacting 
electron Green function in the relativistic case (see Ref. |22] 
p. 32 for Ep = 0) and as well as in the nonrelativistic case 
(up to a factor 7 °, see Ref. |§, p. 124). 

The induced current is —ietr[y G N (x,x)]. But, this 
formal expression is infinite and must be regularized. The 
limit involved in the expression G%(x,x) is discussed in 
detail in section ||: 

G°m(x, x) — — lim lim (G° w (i + z, x — z) 
2 z^o 2 °^o+ v ' 

+G%(x- z,x + z)). 

This definition agrees with experiments in nuclear physics 
[0 and atomic physics (Ref.j2^] p. 275). 
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With this definition, we obtain 



Gj^(x, x) 



^ E n <E F 



E n >E f 



In particular, we obtain the charge density: 
p(x) = J {x) = -ieti[ 1 °G° N {x,x)] 

= |( E i^wi 2 - E i^w 



E n <E F 



E„>E F 



This is indeed the density obtained by Dirac. One 
might wonder how this charge density is related to the 
many-body charge density (or the non-relativistic one), 
which is obtained by summing over the occupied states 
with positive energy. The answer is that we must subtract 
the contribution of the Green function Gq(x, x) (for which 
Ep = 0). Then we obtain 



-ieti[ 1 (G N (x,x)-G° (x,x))] 



* E 

0<E n <E F 



l^n(x)| 5 



The right hand side is the familiar expression for p(x) . It is 
not strictly correct, since it neglects the effect of Gq(x, x), 
which describes the polarisation of the vacuum by the ex- 
ternal potential (i.e. the Uehling potential, see Ref. M| p. 
327 or Ref. [|2| p. 265). It might seem strange that the QED 
and the many-body densities are different, since they are 
given by the same equation 



/oo 
da,tr[ 7 °G&(x,x; 



>)]■ 



The reason for this discrepancy is that the expression for 
p(x) is divergent, and the difference comes from a dif- 
ferent renormalization of it. In the many-body case, an 
innocent-looking convergence factor e lujri is introduced in 
the integral over u (see Ref.j3(| p. 124), which selects pre- 
cisely the states E n < Ep (the states E n < are discarded 
from the start). On the other hand, in the QED case, the 
integration contour is below the real axis from — oo to Ep 
and above it from Ep to oo. 

Several approximations are possible (no-sea, no-pair). 
They are discussed in Refs.§7|, §1. 



5.3 Subsequent works 

This idea of a Dirac sea in an external potential has re- 
ceived many experimental and theoretical confirmations. 
The first and strongest one being the discovery of the 
positron. 




were used as a 
§,@,@,for 
and for rclativis- 



Since then, the Dirac sea and Eq 
basis for strong field electrodynamics 
thermal field theory g§, @> @ 
tic quantum many-body theory in Riemannian spacetimc 

The concept of a Dirac sea is sometimes considered to 
be rather out of date. However, according to Jackiw |j47| , 
it is still the best physical picture to understand striking 



phenomena such as chiral anomalies |48[ and fractional 
fermion numbers ff9|| . 

The case when the Fermi energy is degenerate will not 
be considered in this paper. 



6 The induced current 

It is important to derive the proper expression of the four- 
current as a function of the electron Green function. As a 
field operator, the current must be self-adjoint (it is mea- 
surable), conserved, gauge invariant and should change 
sign under charge conjugation. 



6.1 Definition 

We follow Refs.Jlf 
operator 



|51(1 and define the bilocal field 



2 ^ 

a/3 



When this operator is evaluated at the ground state 
of the system, its value can be expressed in terms of the 
electron Green function. 

The Stueckelberg-Feynman Green functions that we 
use in the present paper are defined from the time-ordering 
operator (@ p.150) 

Ttp a (x)tp (y) = 9(x° - y Q )^ a {x)^ p {y) 
-e{y -x°)j,p{y)ip a {x) 



by 



(Tip a (x)ip (y)) = ihcS a0 (x,y). 



(23) 



We consider here the interacting fields and not the free 
fields as in Ref.(H)) p.272. Therefore, S(x,y) is the full 
electron Green function. 

To define the current we follow Schwinger's prescrip- 
tion ]50| . Let z be a four- vector with z° > 0, the four- 
current is 

J»(x) = -iehcti^Six, x)], (24) 
where S(x,x) is defined by the symmetric limit 



S(x,x) = -\im lim (S(x + z, x — z) 

2 z^O 2 o^ + v 

+S(x — z, x + z)) . 



(25) 



More precisely, we start from a non-zero 3-vector z (i.e. 
|z| 7^ 0) and a positive time z°. We first make z° — > and 
then |z| — > 0. Because of the symmetric limit we obtain 
the same result with z° < (just take z — > — z in the 
above definition) 

S(x, x) = — lim lim (S(x + z, x — z) 

2 z^0 z o^ - v ' 

+S(x — z, x + z)j . 
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From definition (|22|), for z° > we have Moreover, if x and y are space- like separated, we can 

rewrite the bilocal current in terms of field commutators: 

S a0 (x + z,x-z) = ~(ilt a (x + z)j> l3 {x-z)), (26) = J»{x,y)+J»(y,x) 



When we compare this with the definition of the bilocal 
current field operator, we obtain 



4 

a/3 



(29) 



The fact that the current J»(x) defined by Eq.(|24|) 
J (x) = Inn ^ lim + (J (a; + z, x — z)). j s rea j can a i g0 De obtained directly from a property of 

S(x,y). If z is a space-like vector, causality ensures that 



On the other hand, if z° < 



(ip a (x + z)tp^(x - z)f = ipp(x ~ z)i/jI(x + z) 
J^ix) = lim lim (J^tx — z.x + z)). ,\ , \i i \ 



Thus, when the definition of J»(x) is written in terms of Now Eq.(g§) gives 
the bilocal operator, the first argument of the bilocal oper- 
ator is always later than x. It is not clear that the current (ipa(x + z)ip^{x - z)) = ihc^S a \(x + z,x- 2)7°^, 
defined by Eq. (|24|) is real, as a measurable quantity should a 

Now we prove that the current is real and transforms an< ^ §i ves 
properly under charge conjugation. t _ 



6.2 Self-adjointness Thus, we obtain, for space-like z 



Using the definition xp — t/^7 , the bilocal operator can 



S'(x + z,x — z) = — 7 S(x- z,x + z)j°. (30) 



be rewritten Here the dagger operator acts only on the spin- variables: 

a p(x+z,x-z) = Sp a (x+z,x—z). From Eq. (^p|) and the 



J>(x, y) = - ty(x)h°rm ~ ^XtVJ^W] , identity 7 o 7 , 7 o = y 4 we can check that j, {x y = j, [x) _ 
where A T is the transpose of A. From the identity (7°7 Al )^ = 

7°7 A ' (Ref . fiof p. 693) we obtain immediately 6.3 Charge conjugation 

J M (x, yY = J M (y, x). (28) Let C be the charge conjugation operator acting on fields 

and C = i7 2 7° the charge co njug ation matrix. We have 

Now we can make an important remark concerning the the identities ((To) p.l52g and j53| p. 70) 
bilocal current. If x and y are separated by a space- like 

interval (i.e. if (x — y) ■ (x — y) < 0), then ip(x) and i/>(y) Cip(x; e)C — Cip(x; — e), 
anticommute. This property is called causality and means Cxp(x; e)C^ — tp(x; —e)C. 
that two events that are too far to be linked by a light- 
ray are independent. Causality is a basis of all axiomatic In the presence of an external field, the sign of the charge 
approaches to quantum field theory Jjj^]. in the field operators is reversed under charge conjugation 

Thus, if (x - y) ■ {x - y) < we have (H P-19)- 

The action of the charge operator on the bilocal cur- 



J M (^ = + i> a {y)M x )) rent is 

CJ»(x,y;e)tf = ^( x; - e )C^C^(y;-e) 



-rP(x ] ~e)(CrC) T xP(y;~e)). 

Combined with Eq.(p8|), this means that the bilocal cur- 
rent is self-adjoint for space-like separated points. From the identities (Q p.693) C^C^ = -7^ and C T = 

The symmetric definition (25) of the current in terms C = C = — C we obtain 

of the Green function was chosen such that first z° — > + . „. . f . 

with |z| ^ 0. This is to ensure that, for z° small enough CJ K x >y> e ) c = J \x,y;-e), 

(i.e. z° < |z|), the interval Z becomes space-like and the 2 Their misprinted definition of C is corrected on page 693 

current is real. of lEql. 
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Therefore, the sign of the bilocal current is reversed under 
charge conjugation. 

Now we investigate the behaviour of the current un- 
der charge conjugation, when it is defined in terms of the 
Green function. To do this, we need an equation for the 
Green function of the charge conjugated problem. We ob- 
tain it by transposing both sides of Eq.(17]). This trans- 
position is meant for the 4x4 matrices only: S T ss '(x, y) = 
S s ' s (x, y). Then we multiply both sides of the equation by 
the charge conjugation matrix C on the left and C _1 on 
the right. We use the property C^ T = — 7 M and we 
obtain 

5{x — y) = (ihcj ■ d — rac 2 + ej ■ A(x))S c (x, y) 
SS c (x,y) 



where we have defined S c (x,y) = CS(y,x) T C~ 1 . If we 
compare this equation with the second basic equation ( |l3| ) , 
we see that they become identical if we make the trans- 
formation e — > — e, or equivalently — >■ — and j M — ► 



— j M . But, using Eq.fllOD and the charge conjugation prop- 
erty of the induced current J' i (x), we see that j M — > —j 1 * 
implies A^ — > —A^. In other words, we have the symmetry 

S c (x,y) = S(x,y--j) = CS(y,x;j) T C-\ 

Therefore 

tr[^S{x ~z,x + z; j)} = tr[Cj fl C- 1 S(x + z,x-z; -j) T ] 
= —tr[-f^ T S (x + z,x — z; — j) T ] 
= — ti[j^S(x + z,x — z; — j)]. 

And the current can be rewritten as the antisymmetric 
limit 

S(x, x; j) = — lim lim (S(x + z.x — z: j) 

-S(x + z,x- z;-j)). 

This equation shows that the current is an odd function 
of the external current: J M (x;j) = — J^x; — j). Changing 
the sign of j amounts to changing the sign of A, and we 
can also write Jp{x; A) = —Jp(x; —A): the current is an 
odd functions of A. Now, since the electron loops of QED 
are generated by the current, and since the x in S(x, x; A) 
is linked to another photon propagator, we see that the 
electron loops with an odd number of photon lines are 
zero. This is a version of Furry's theorem for QED with 
external field. The relevance of charge conjugation for the 
QED current was already noticed by Kramers in 1937 j54j. 
Related results can be found in Ref.0, sections 4.2, 9.4 
and 15.1. 



6.4 Charge conservation 

Charge conservation is described by the equation d li J^{x) = 
0. How does this translate for our bilocal operator? We 
want a real and symmetric expression for d^J^ix) = 0. 



Why symmetric? If a bilocal operator F(x, y) is antisym- 
metric (i.e. F(x, y) — —F{y, x)), then its value at x — y is 
zero, except at divergences, such as (x — y) / (x — y) 2 . Since 
the regular part of F(x, y) is zero on the diagonal, the 
use of a symmetric operator eliminates the antisymmet- 
ric divergences without changing the value of the regular 
operator on the diagonal. 

An obvious candidate is (9? + df l )3 ,J '(x, y)/2. From the 
field equations p. 195 

(ihcdf, - eA^^ip - mc 2 xp = 0, 
{ihcd^ + eA^xjj^ + mc 2 4> = 0, 

we check that (0£ + 8^)3^ (x, y) = 0. 

In terms of the electron Green function, this becomes 

zcfic 

—dftr[~f' t S(x + z, x - z) + ^S{x + z, x - z)] , 

where we have used df. +z + df~ z = df.. and where z is 
space-like. Charge conservation can then be directly veri- 
fied from Eqs.@ and (0). 



6.5 Mass-reversal symmetry 

The concept of mass-reversal symmetry was discussed by 
Peaslee 55 and Tiomno pqj. For massless particles, mass- 



reversal symmetry is known as chiral symmetry. 

The mass-reversed fermion field i/>(— m) satisfies the 
same equations as the usual fermion field ij)(rn), but with 
a reversed mass: 

(ihcry ■ d — e7 • A + mc )?/>(— m) = 0. 

Since 7^75 = —757^ we can take i/)(x; —to) = 7 5 ^>(x;m) 
(up to a phase factor r\ m which disappears in the Green 
function) . To see how the Green function transforms under 
mass reversal, we multiply both sides of Eq.([l3|) by 75 on 
the left and on the right. The commutation rules between 
75 and 7^ yield 

S(x — y) = (—ihay ■ d — mc 2 + ■ A(x))~/5S(x, y; 771)75 

+l eh n , l5 —y- — 75, (31) 

which is the equation for —S(x,y;—m). Therefore, 

S(x,y;-m) = -j 5 S(x,y;m)j 5 , 

and the mass-reversed current is 

J fJ '(x; — m) = —iehctT[ r y fJ 'S(x 1 x; —to)] 

= ie/ictr[757 M 75S'(a;, x; m)] = J^(x;m). 

Thus, the induced current is even under mass-conjugation, 
and the vector potential also, as expected [ p6| . 

In other words, the induced current is an even function 

of TO. 
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6.6 Gauge invariance 

The Lagrangian of QED is invariant under the gauge trans- 
formations (ficj p. 64) 

tj, A ( x ) = e-^^xPix), 
A A {x)=A li {x)+Kcd ll A{x). 

After a gauge transformation, the bilocal current for space- 
like separated points Eq.(^) becomes 

J A (x, V) = cos{eA(x) - eA(y))J(x, y) 
+ — sin(eyl(a:) - eA(y)) 

In terms of the electron Green function, we must consider 
the transformed Green functions 

S A (x,y) = e~ %eA ^ +leA ^S{x,y). (33) 

Hence, 

S A (x + z,x - z) + S A (x - z,x+ z) = 
cos(eyl(a;) — eA{y))[S{x + z,x — z) + S(x — z, x + z)) 
—i sin(eA(x) — eA(y)) (S(x + z,x — z) — S(x — z, x + z)) . 

A priori, the singularities of the Green function are not 
weak enough to make the second term tend to zero when 
z — > 0. Therefore, our definition of the current is not ob- 
viously gauge invariant. 

This problem can be solved by using a Green function 
multiplied by Dirac's phase [^9| 

/i 
dXz^A^x + Az)] 
-i 

xS(x + z, x — z). 

A generalization of this phase was used recently for the 
renormalization of infrared divergences in QED p?] ] , J58| . 

The modified Green function S'jx + z,x — z) is gauge 
invariant because, according to Eq.(p3) in the appendix, 

J d\d l _ l A(x + Xz)z^ = A(x + z) - A(x - z). 

The right hand side compensates for the gauge transfor- 
mation phase in Eq. (|33|) , so that S' A (x, y) — S'(x, y). 

If we substitute S'(x + z, x— z) for S(x + z, x — z) in the 
definition of the current (|2J), we obtain a current which 
is real, odd under charge conjugation and gauge invariant, 
but charge conservation becomes doubtful. Thus, we stick 
to the original definition Eq.(|24|). 

7 Renormalization 

The renormalization of QED in the presence of an external 
field does not seem to have been treated in detail. Several 



attempts exist jjj, |5j|, [[59], pl| , but none of them was 
fully completed. We start now a renormalization of QED 
in an external field which closely follows the renormal- 
ization of vacuum QED. The first step is to determine 
Dyson's relations, which are non perturbative expressions 
that link rcnormalizcd and non renormalized propagators 
and fields. Then we give a non perturbative relation be- 
tween the renormalized potential and the renormalized 
vacuum polarization. Finally, we describe the renormal- 
ization rules that enable us to obtain finite quantities for 
any strong field or SCF Feynman diagram. 

7.1 Dyson relations for strong field QED 

The first point is to determine the Dyson relations between 
renormalized and unrenormalized propagators. This rela- 
tion is given by Sterman |6(i[| for scalar fields and we derive 
it for QED. If we denote the renormalized quantities with 
a overlined symbol, in the presence of an external poten- 
tial a^(x), Dyson's relations are |HJ 

S(x,y;a;e ,m ) = Z 2 S(x,y; \/z^a;e,m), 
D(x, y; a; e , m ) = Z 3 D(x, y; \fZ3~a; e, m), 
A(x; a; eo, mo) = \fz~ 3 A(x\ \fz~ 3 a\ e, m). 

To derive these relations, we expand S over a: 

b[x,y;a) = S(x,y;0) + J dzi — j—^ a^zt) H 

The Dyson relation for the vacuum propagator S(x, y; 0) 
is S(x, y; 0) = Z 2 S(x, y; 0). For the other terms, the func- 
tional derivative with respect to the external potential 
brings a factor A(x) in the path integral. Therefore, 

5 n S{x lV -0) _ r7n/2 S n S(x,y;0) 

— ^2^1 ; ■ 

Sa n 6 5a n 

Hence, 

S(x,y;a) = Z 2 S(x,y;0) + Z 2 [ dzi ^^' y ' °^ y/z^a^zx) 
H 

= Z 2 S(x,y; \[z~ 3 a). 

Or, more precisely, 

S(x, y; a; e , m ) = Z 2 S(x 1 y; \fz~ 3 a\ e, m). 

Because of this relation, we also define the renormal- 
ized external potential as a M (x) = ^fZ^a^ix). It may seem 
strange that the external potential (x) and the full po- 
tential A^x) are not renormalized with the same formula. 
However, we must recall that a^x) was defined in terms 
of an external current j^x), which was used to generate 
Green functions by functional derivatives. The bare Green 
functions are generated by A^ (x) (x) and the renormal- 
ized ones by A^ (x)j M (x) . Since we want these two terms to 
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be equal, the relation yfZ^A p (x) = A p (x) implies J Ax) — 

VZ^j^x) (see Refs.[|o) p. 297 and |l9) p.288). More phys- 
ically, we can say that an external field is made by prepar- 
ing a certain density of (charged) matter n(x). Before 
renormalization, this density corresponds to a charge den- 
sity j°(x) = erjn(x). The renormalization modifies the 
charge, but not the density of matter. Thus the renor- 
malized charge density is now j°(x) = en(x) = y/Z~3j°(x), 
and the corresponding relation for the external potential 
is again a M (x) = \fZ^,a^(x). 

In terms of the renormalized propagators, the func- 
tional equation becomes (using eo = e/^/Z^ and mo = 
to + 5m) 

(iHcy ■ d — mc 2 — Smc 2 — ej ■ A(x))Z2S(x, y) = 5(x — y) 



-iehcZ2j' J ' / ds 



SS(x, y) 
8a\(s) 



D%(s,x). 



and 



ZzA^x) = a p (x)/ V^3 - i—f=hcZ 2 

V Z3 

x / d S DUx,y)tT[YS(y,y)}, 



or 



A, l ( x ) = ^.-i^hc I ,lsD' l ] l A.r. !l nv[~."S [! ,. ! ,)}. 



(34) 



7.2 Dyson relations for SCF-QED 



Notice that, in the case of SCF-QED, Dyson's relations 
are much simpler, because e$A = eA: 

Six, y; e A; e , m ) = Z 2 S{x, y; eA; e, to), (35) 
D(x, y; e A; e , m ) = Z 3 D(x, y; eA; e, to), (36) 
A(x; eoA; eo,mo) = \fz~lA(x;eA;e,m). (37) 

The main difference between vacuum QED and SCF- 
QED is the renormalization of the induced vector poten- 
tial. The total bare potential is 

A p (x) = a p (x) - ie hc J dsD° v (x, s)tr[y v S(s, s)],(38) 

and the problem is the renormalization of tr[7 l/ S'(s, s)]. To 
indicate that the electron Green function is calculated in 
the presence of the full potential A p (x) we denote it by 
S(x,y;A). 

Dyson's relation ( p5|) gives us 

tY[YS(s,s;A)} = Z 2 tT[ 1 v S(s,s;A)}. 

We rewrite the right hand side using Eq.(]5l|): 



tr[-fS(s, s;A)}=tT[j v S(s,s;0)] 
u 5S(s, S ;\A)- 



dX / dx 



xtr 



Now — iehctv[ r y v S{s, s; 0)] is the current in the presence 
of a zero external field A = 0. Thus, it is the vacuum cur- 
rent, which is zero (more precisely, trpy" S^s, s; 0)] is made 
of fermion loops with one external photon line, which are 
zero by Furry's theorem). 

The second term is also well known from the renor- 
malization of the vacuum polarization (see e.g. Eq.(22) of 
Ref.flel) 



iehcZiix 



5S{s, s;XA) 



5A»(y) 

n^(s,y;XA) + (Z 3 - l)D°^( S ,y). 
Thus, we obtain 



(39) 



iehcZi J dytr{jvS(s, s; A)] = 



dX / dyn u Js,y;\A)A»(y) 



+(Z S -1) J dyD\Us,y)A^y). 

After this reorganization of the induced current, Eq.([38|) 
becomes 

A p (x) = a^(x) + 3 _ A^x) 



+ - 



dsdydXD°(x, s)n»»(s, y; XA)A p (y). 



Now we know that the relation between the bare and 
renormalized potentials is A p (x) = \[ZzA p {x). This yields 



A p (x) = a^x) + J dXj dsdyD%(x,s) 
xn»P( S ,y;XA)A p {y), 



(40) 



where the true external potential a p (x) is renormalized by 
a p (x) = y/Z^a^x). 

It should be noticed that, in all our manipulations, 
the renormalization factors Z2, Z% and 5m where taken 
at A^ — 0. In other words, the renormalization factors of 
QED with external field are the same as the renormaliza- 
tion factors of vacuum QED. More precisely, the vacuum 
renormalization factors removes the divergences of QED 
in external field. However, renormalization conditions may 
introduce finite differences between the renormalization 
factors of vacuum QED and QED in an external field. 

From Eq. (^0|), it can be checked that 77 is indeed the 
renormalized vacuum polarization. If we write X p {x) for 
the second term on the right hand side of Eq. ([lo]) we find 



D^{x,y) = 



5A p (x) _ 5a p (x) SX^x) 

ST(y) ~ 5T(y) 5-f(y) 

6X p (x) 5A p {x) 
Z 5A p (z) 5T(y) 

<U / dzD° pa (x, S ) 



= DUx,y) + 



xn°t>( S ,z)D p „(z,y), 



(41) 
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dsD%(x,s)n°»( S ,z). 



where we used 

6AJz) 



This last equation is proved by using Eq.(|39|) to rewrite 
Xfj,(x) in terms of a functional derivative of S, and by 
using Eq.(^J) to carry out the integral over A. 

Equation (^l|) expresses the usual relation between the 
photon Green function and the vacuum polarization. 



8 Weak field renormalization 

In this section, we give the renormalization rules of vac- 
uum QED, which are valid for weak field renormalization. 
We repeat some of the rules given in for complete- 
ness, and because we work now in the direct space (and 
not in the Fourier space). The first set of rules is equivalent 
to Zimmermann's forest formula for the removal of subdi- 
vergences. The second set of rules subtracts the superficial 
divergences. 



8.1 Removal of subdivergences 

If r is a Feynman diagram, and U(r) is the regularized 
value of the diagram as a function of external momenta 
and masses and regularization parameters, then the value 
of the diagram with subdivergences subtracted is R(r) 
defined by @, @ 

{71,7V,.. . > 

xc/or/h^^,...}). (42) 

In this expression, the sum runs over all sets of disjoint 
renormalization parts of r. A renormalization part 7 of r 
is a one-particle irreducible (1PI) subgraph of r, different 
from r itself, such that 7 has two or three (amputated) 
external lines. A diagram is 1PI when it is connected and 
cannot be disconnected by cutting through any of its inter- 
nal lines. Two renormalization parts 7 and 7' are disjoint 
if they have no vertex in common. 

There are three types of renormalization part: self- 
energy 1PI diagrams (i.e. 1PI diagrams with two am- 
puted external electron lines), vacuum polarization 1PI 
diagrams (i.e. 1PI diagrams with two amputed external 
photon lines) and reduced vertex 1PI diagrams (i.e. 1PI 
diagrams with two amputed external electron lines and 
one amputed external photon line). The index (i) depends 
on the type of the renormalization part. If 7 is a self energy 
1PI diagram then i is or 2, if 7 is a vacuum polariza- 
tion 1PI diagram then i is 3, if 7 is a reduced vertex 1PI 
diagram then i is 1. 

Finally let us define T '/ {7(1), lu/y ■ ■ ■}■ We start by the 
definition of r/j^ . It varies with the type of renormal- 
ization part and with the index (i). 

If 7 is a self-energy 1PI diagram, then i can be or 
2. If the complete diagram is r and 7 a self-energy 1PI 



subdiagram of r, then in the term 1^/7(0) , the 1PI diagram 
7 is replaced by a point and in the term r/j^) the 1PI 
diagram 7 is replaced by a free electron line. 

If 7 is a vacuum polarization 1PI diagram, then i is 3. If 
the complete diagram is r and 7 is a vacuum polarization 
1PI subdiagram of JT, then in the term r/jrss, the 1PI 
diagram 7 is replaced by a free photon propagator. 

If 7 is a reduced vertex 1PI diagram, then i is 1. If 
the complete diagram is r and 7 is a reduced vertex 
1PI diagram, then in the term r/jn\, the vertex dia- 
gram is replaced by a point. The terms -T/{7(j) , 7(jn , . . . } 
are then defined recursively. For instance, to define r" = 
r /{7(i),7('i')}' we nr st put F' = r/rf(i), so that T" = 

We insist here on the fact that the removal of subdi- 



vergences given by Eq.(42) is valid for any (connected or 
disconnected) Feynman diagram. It is the expression, for 
each diagram, of Dyson's relations. The reason why the 
counterterms 6*1(7) are required, although Z 2 = Z\, is 
explained in fl62| . 

In Eq.(p2[) all counterterms Ci(j) are assumed to be 
known from the renormalization of the superficial diver- 
gence of the subdiagram 7. 



8.2 Superficial divergences 

In the second step, we determine the counterterms C'i(r) 
of the divergent graph r. If the diagram r is not 1PI, 
the counterterms Cj(-T) are zero. If r is 1PI, we must 
distinguish three cases. A self-energy diagram is linearly 
divergent, thus we must remove two terms. From Lorentz 
covariance, we can write the renormalized value of the 
diagram r as 



R{T; x, y) = R(r; x, y) + C (r)c 2 S(x - y) 

+c 2 (r)s°-\x, y ). 



(43) 



A vacuum polarization diagram is quadratically diver- 
gent, thus we should have to remove three terms. However, 
Lorentz covariance and the Ward identities for the pho- 
ton propagators cancel the first two counterterms, and the 
renormalized value of r is 

R^{r-x,y) = R^{r-x,y) + C z {r)D Q -l{x lV ). (44) 

Finally, a reduced vertex diagram is logarithmically diver- 
gent, and its renormalized value is 

R(r; x, y; A, z) = R{r- x, y: A, z) 

+C 1 (r) 1 x 5(x-y)8(y-z). (45) 

The infinite constants C, (T) are determined from the renor- 
malization conditions. 

It was noticed in Ref.|62| that 



c 2 (r) 



(46) 
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where Fj is the vertex diagram obtained by branching a 
photon line on the j-th free electron propagator of the 
self-energy diagram T . 

Since the renormalization of weak field QED is exactly 
the same as the renormalization of vacuum QED, we refer 
the reader to Refs.|lC| and j6^] for examples. 



9 Strong field renormalization 

For the strong field renormalization, the removal of sub- 
divergences is again given by Eq.(|42j), which is a purely 
algebraic identity related to the Hopf algebra structure of 
renormalization Q, p5[ . 

According to the general strategy of renormalization, 
the non-renormalized diagrams are evaluated in terms of 
renormalized electron propagators. In vacuum QED, the 
mass m used in the free electron propagator is finite, it is 
not the infinite mass mo of the bare free electron propa- 
gator. Similarly, the charge and the potential used in the 
basic electron propagator of strong field QED are finite. 
In other words we take 



Using Eq.(45), the counterterm for the second diagram is 



S a = (ihcj ■ d 



rac 



■a(x))~ 



Therefore, the external potential diagram means now 

U(- — *) = a M (x). 

To remove the superficial divergences we start by ex- 
amining all possible superficially divergent diagrams (Chap- 
ter 8 of Ref . |l^] ) . A priori, the superficially divergent 1PI 
diagrams are the self-energy and vertex diagrams, and the 
1PI diagrams with one external photon (tadpole), two ex- 
ternal photon (vacuum polarization), three external pho- 
tons and four external photons (scattering of light by 
light). 



9.1 Self-energy and vertex diagrams 

To renormalize self-energy diagrams, we take the exam- 
ple of the one-loop diagram, and we expand the electron 
propagator using the Born expansion (fl6q| p. 572) 



S a = 5° + eS°a • 7 S l 



+ e 2 S°a ■ 



~fS°a-~fS a . 



Thus 



The strong field diagram on the left hand side is denoted 
r. The first diagram on the right hand side (denoted r') is 
a self-energy diagram of vacuum QED, the second diagram 
(denoted r{) is a vertex diagram of vacuum QED and the 
third diagram is finite by power counting. Using Eq. ([hf) , 
the counterterms for r are 

R{r'; x, y) = R(r'; x, y) + C a (r')c 2 S{x - y) 

+c 2 (r')s n -\x iy ), 



dzR(r[; x, y; A, z)a\{z) = J dzR(r{; x,y; A, z)a\(z) 

+C 1 (r[) 1 x S(x-y)a x (y). 

The total counterterms for the strong field diagram r are 
obtained by adding the counterterms for i~" and r[. From 
the relation ( |46| ) between C\(r[) and C^iT') we obtain 
for the strong field self-energy diagrams 

R{r- x, y) = R{r- x, y) + C (r)c 2 5(x - y) 

+C 2 (r)(S°~ 1 (x, y)-ery- a(y)S(x - y)) 
= R(r;x,y) + C (r)c 2 6(x~y) 

(47) 



-C 2 (r)S a \x,y), 



where C (r) = C (r') and C 2 (r) = C 2 {r'). Therefore, 
the strong field self-energy diagrams are renormalized with 
the same formula as the vacuum diagrams, the countert- 
erms Co(-T) and Ci(-T) are the same as for the vacuum 
case, the only difference is that the free propagator S° is 
replaced by the strong field propagator S a . More precisely, 
the counterterms are the same as those of vacuum QED if 
the renormalization conditions are compatible. Otherwise, 
their difference is finite. 

For a general strong field self-energy diagram r, the 
proof is similar. The subdivergences are removed with for- 
mula ([42]). The same Born expansion is made. The first 
diagram r' is now the diagram _T (with subtracted sub- 
divergences) where all strong field electron propagators 
S a are replaced by free electron propagators S , and 
is the diagram obtained by adding a photon line to the 
j-th electron line of r'. The remaining diagrams are fi- 
nite by power counting. The superficial divergences of r' 
and -Tj are removed by the vacuum QED renormalization 
prescription, and we find Eq.(|47|) again. 

For a strong field vertex diagram, we have a logarith- 
mic divergence and the result is exactly the same as for 
vacuum QED: 



R(r-,x,y;X,z) 



R(r-,x,y;X,z) 

+d(r) 7 x 6{x- 



y)5(y-z), (48) 



where Ci(P) is the same as for vacuum QED, up to pos- 
sible finite terms. 



9.2 Vacuum polarization and tadpoles 

We still have the subdivergence-subtracted 1PI diagrams 
with one, two, three or four external photon lines (and no 
external electron line). If the diagram has three or four 
external photon lines, then it is finite. To prove this, we 
expand the strong field electron propagator using the Born 
expansion. Then we are back to the case of vacuum QED, 
where we know that a 1PI diagram with three external 
photon lines is zero and a 1PI with four external photon 
lines is finite (To) (once subdivergences are subracted). A 
1PI diagram with two external photon lines is a vacuum 



16 



Christian Brouder: Renormalization of QED in an external field 



polarization diagram, which is renormalized as for vacuum 
QED with Eq.@ 

R^(r; Xl y) =R fiV (r-x,y) + C 3 {r)D°-l(x,y), (49) 

where Cs(r) has the same value as for the corresponding 
vacuum QED diagram, up to a possible finite quantity. 

Now comes the diagram which is really typical of QED 
with external field: the tadpole. The induced vector po- 



tential A^(x) 



L {x) is given by the sum of all tadpole 



diagrams r. In vacuum QED, the induced vector poten- 
tial is zero and the tadpoles do not intervene. For the case 
of QED with an external field, once subdivergences are 
removed, the tadpole has a cubic superficial divergence. 
But because of Furry's theorem (i.e. of the fact that the 
current is odd under charge conjugation), the only diver- 
gence that remains is the same as for vacuum polarization. 
We illustrate the general procedure with the simplest ex- 
ample (see also [|2j). By using the Born expansion for S a 
we obtain 





Furry's theorem eliminates electron loops with one and 
three external photon lines (the first and the third dia- 
grams on the right hand side). More precisely, since the 
induced current is an odd function of the external poten- 
tial, the first and the third diagrams are not present on 
the right hand side. The last diagram is finite, because 
if we make a further Born expansion of S a in it, we ob- 
tain a vacuum photon-photon scattering diagram, which 
is finite, plus a diagram with five external photon lines, 
which is finite by power counting. Therefore, the only di- 
vergent diagram on the rhs is the second one, which is a 
vacuum polarization diagram. We denote r{ this vacuum 
polarization diagram and r the strong field tadpole on the 
lhs. 

We can write the bare potential as a sum of Feynman 
diagrams 

A x (x) = a x (x) + Y,e l o rl Ux(x;r), 
r 

where \r\ is the number of vertices of the tadpole diagram 
r. If r is the simple tadpole of the above example, our 
discussion shows that its value is 

U x (x;r) = e J dy&zDl^x^R^ir^-y.z^iz) 

+finite terms. 

Now, we introduce the counterterm of r{ using Eq.([4~4|), 
and we obtain the renormalized induced potential 



Rx(x;r) = e J dydzDl^x^R^iri-y^^z) 
-\-eCz(r[)a\(x) + finite terms. 



Similar equations can be found in j66[ p. 552, |5Sj ] . 

Therefore, the counterterm for the tadpole diagram r 

is 

Rx{r-x) = R x {r-x) + eC 4 (r)a x (x), 

where C 4 (r) = C 3 (r{). 

For a general tadpole diagram r, the result is similar. 
If Rv(r;x) is the value of r when all subdivergences are 
removed, the superficial divergence is removed with the 
counterterm 



where 



R»{r- x) = R*{T\ x) + eC 4 {r)a^{x), 



c 4 (r) = ^c 3 (r i ). 



(50) 



In the last relation, 1^ is the vacuum polarization obtained 
from the tadpole r by adding a photon line to the i-th 
electron line of r and by transforming all electron lines 
S a into free electron lines S° in the resulting diagram. In 
other words, Fi is a vacuum QED vacuum polarization 
diagram, and C^-Tj) the corresponding counterterm. 

After this analysis, we obtain the following addition 
to the description of the renormalization of vacuum QED. 
There is now a new index i = 4 and new IPI diagrams 7 
which are tadpoles. In the term 7^/7(4), the tadpole sub- 
diagram 7(4) becomes an external potential line »~~*. 



9.3 Summary 

In this section, we summarize the renormalization rules for 
strong field QED. For any strong field Feynman diagram 
r, the subdivergences are subtracted by applying the for- 
est formula ( 42]) . The only difference with the weak field 
case is that the vacuum electron propagators are replaced 
by strong field electron propagators. 

When the subdivergences are subtracted, the superfi- 
cial divergence is treated as follows. If r is not IPI, the 
superficial counterterms are zero. In particular, if r con- 
tains tadpoles, then the counterterms Ci(-T) are zero. This 
is because a tadpole inside r can always be disconnected 
from r by cutting a photon line. 

If r is IPI, then the superficial divergence is removed 
using Eq.p7|) for a self-energy and the indices are i = 
and i = 2, Eq.(48) for a vertex and the index is i = 1, 
Eq.(^9|) for a vacuum polarization with the index i = 3 
and Eq.(po|) for a tadpole with the index i = 4. 



Now we give a few example to illustrate the general 
rules. 



9.4 Example 1: the potential 

As far as we know, the renormalization rules for strong 
field QED are new. Therefore, a number of examples are 
required to see how they work in practice. 
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We start with the renorm alizat ion of the potential for 
the diagrams given in section 4.2.1. The first diagram has 
no subdivergence. It is renormalized with Eq.(|50|). 

The second diagram is not 1PI. To renormalize it, we 
renormalize its two disjoint divergent subdiagrams: the 
vacuum polarization diagram on the left (index i = 3) 
and the tadpole on the right (index i = 4), and we multi- 
ply the renormalized diagrams. The general rule for a non 
1PI diagram _T is the following: write J 1 as a product of 
1PI diagrams, renormalized each 1PI diagram, and multi- 
ply the renormalized diagrams Jl8| . This rule is equivalent 
to the forest formula (E2h for non 1PI diagrams, because 
renormalization parts belonging to different 1PI subdia- 
grams are disjoint. 



R 



= U 



+e 



For the last diagram, there are two subdivergences: on 
the left, a self-energy 7 (indices i = and i — 2) and on 
the right, a vertex 7' (index i = 1). Since these subdia- 
grams are not disjoint we obtain, adding the superficial 
divergence counterterm 



R 



= U 



+ eC 4 



+C (^)[/(~~ 0) + C a (£*)U{~- Q) 

+Cl (^)u(~~Q). 



To check these results, we make use of the renormal- 
ization factors. 



Zi,n is also the sum of all — C^-T) where r runs over the 
(1PI) self-energy diagrams with n loops, Z 3 .„ is the sum 
of all Ca(r) where F runs over the (1PI) vacuum polar- 
ization diagrams with n loops, Z 3-n is also the sum of all 
Ci(r) where r runs over the (1PI) tadpole diagrams with 
n loops. Finally, 5m n + Ylk=i $ m kZ2,n-k is the sum of all 
Co(-T) where r runs over the (1PI) self-energy diagrams 
with n loops. 

For our example, we have 

Thus, the renormalized induced potential up to two 
loops is 



eU 



e 3 U 



z 3 U 



-(e 2 Z 3 ,i + e 4 (Z 3 , 2 + Zl !))C/(— ^ 
-e 3 Z 3 ,it/(— O) +e 3 Z 3 , 1 C/(^0^) 



-e 3 5m 1 U 



To check this result, we take Dyson's relation between 
renormalized and bare electron Green functions 

Z2(e, m)S(x, y; ea; e, m) = S(x, y; — ; , m + 8m). 

Z3 \JZ-i 

We expand the right hand side 

Z 2 (e, m)S(x, y; ea; e, m) = E/(=) + e 2 Z 3il U(^*) 



-e SmiUi 



e 2 U( 



Zo = 



1 + J2 e2nZ 2,n 
n=l 
00 

1 - ]T e 2n Z; 

n=l 

00 

5m = e 2n 5m 



Z. 



3,n ) 



n=l 



The relation between the renormalization factors and the 
counterterms is as follows |62j . Z2, n is the sum of all C\ (r) 
where r runs over the (lPlJvertex diagrams with n loops, 



+e 2 U( 



where the electron propagator = is calculated in the field 

of the renormalized external potentials a(x). If we use this 
expansion in Eq. (|34| ) we obtain the renormalized potential 
up to two loops. We recall that 



dzS a (x,z)S a (z,y) 



dS a (x,y) 
dm 
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Obviously, it is much simpler to renormalize directly 
from the Dyson relation. However, this does not renormal- 
ized all Feynman diagrams. Only the sum of them is finite, 



as was discussed in |62|. In practice, all the Feynman dia- 
grams of a certain order cannot always be calculated, and 
it is necessary to renormalize each diagram separately. 



9.5 Example 2: the self-energy 

For QED in an external field, the bare and renormalized 
self-energies are defined by 

S^ 1 = ihcj ■ d — mo — eo7 • A — S, 
S^ 1 = iHay ■ d — m — ■ A — S. 

From Dyson's relations between S and S we obtain the 
following relation between £ and S. 

S(ea; e, m) = Z2S{ea/Z^] e/yZa, m + 8m) + Zi8m 
+ (Z 2 - l)(-iHay -d + m + ej- A). 

In the self-energy, the only new diagram with respect 
to vacuum QED is 




This diagram is not 1PI. In vacuum QED, all self-energy 
diagrams are 1PI, but this is not the case in the presence 
of an external field. This diagram is easily renormalized 
as a product of two 1PI diagrams: 

eC 4 (~©)[/(<9l) 

+<?i(<^)E/(~©) +eC 1 («££)C 4 (~©)E/(~^. 



10 Self-consistent field renormalization 

The case of self-consistent QED is simpler, since we have 
the same diagrams as for vacuum QED (with S° replaced 
by S A ). Therefore, the renormalization rules are the same, 
except that the free fermion lines become the fermion lines 
in the SCF potential A(x). 

Since there is no tadpole in the fermion propagator of 
self-consistent QED, the only point which is delicate is the 
renormalization of tadpoles to calculate the self-consistent 
potential. When all subdivergences are removed, the su- 
perficial counterterm is obtained by expanding the SCF 
potential in terms of strong field diagrams and by renor- 
malizing those. The result is 

R"(x; T) = fl"f>; r) + eC 4 (r)A"(x). 

Thus, the counterterm is proportional to the renormalized 
potential itself. 



11 Renormalization conditions 



The renormalization theory is an unambiguous method to 
remove the subdivergences of a Feynman diagram. How- 
ever, once all the subdivergences are removed, we must 
still specify the value of the superficial divergence. This 
is done by giving renormalization conditions. In the case 
of QED without external field, these conditions are well- 
known ( fTcjl p. 413). However, for QED with external field, 
they have never been stated precisely. It is the lack of 
proper renormalization conditions that yields the ambi- 
guities in the results of Dosch and Muller J5lj . There are 
three kinds of renormalization conditions: for the current, 
for the photon propagator and for the electron propaga- 
tor. Some of them were investigated in the early days of 
QED pi. 



11.1 The current 

The renormalization condition for the current is deduced 
from the neutrality of matter. The total charge of a piece 
of matter is the sum of the proton charges and the elec- 
tron charges. Thus, the induced charge density due to the 
vacuum polarization integrates to zero. If this were not 
true, a piece of matter with an equal number of electrons 
and protons would have a net charge. 

In the case of strong field QED, we operate the renor- 
malization rule for tadpoles (| 
obtain 



with D° (x,y) and we 



dyD°"(x,y)R(r;y,v) 



&yD\l{x,y)R{r;y,v) 



+eC A (r)-Ux). 



The left hand side represents the renormalized current in- 
duced by the external current (x) . Two points must be 
noticed here. Firstly, the superficial divergence of the cur- 
rent diagrams (tadpoles) are removed by a counterterm 
proportional to the external current, and any induced cur- 
rent proportional to the external current cannot be dis- 
tinguished from a renormalization of the charge (see [ |66| 
p. 553). With some hindsight, this point can be recognized 
in Ref . |67|] . Secondly, the integral of the induced current 
must be zero to ensure matter neutrality. Thus 

JdxJ dyD°~l(x,y)R(r;y,p) = 

This second point was clearly made by Uehling j33| . When 
the external current is not neutral, this is enough to de- 
termine c^r). 

In the case of self-consistent QED, the counterterm is 
proportional to A(x), and for a neutral atom the corre- 
sponding J{x) integrates to zero. Therefore, the criterium 
of matter neutrality is not enough to determine the renor- 
malized current of self-consistent QED. 
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11.2 The photon Green function 

To specify a renormalization condition for the photon Green 
function, we use the requirement that, very far from the 
external field, low frequency Compton scattering should 
be given by the Thomson formula p(| p. 640 and two 
close charges should interact with a Coulomb potential 
@ p.670 and |o) p.325. 

In other words, for very large x and y and small x — y, 
D^ v (x,y) should tend to D^ v (x,y). 

11.3 The electron Green function 

For the electron Green function, the renormalization con- 
ditions do not seem to have been studied beyond vacuum 
QED. One could probably assume that, for very large x 
and y and small x — y, S(x,y) should tend to S A (x,y). 
We intend to study the validity of these renormalization 
conditions in the specific case of atomic physics. 



12 Conclusion 

In this paper, the equations of self-consistent quantum 
electrodynamics were derived from the Schwinger approach. 
These equations are valid for any number of electrons in 
the system and the functional derivative of the Green func- 
tion in the presence of bound electrons was used to calcu- 
late various perturbative solutions of the Schwinger equa- 
tions. 

An interesting application of self-consistent QED is to 
see if the solitonic solutions of the Maxwell-Dirac equa- 
tions p§| ] survive when the two loop interaction is turned 
on. In other words, can an interacting electron be bound 
by its own vacuum polarization ? 

The general renormalization rules for strong field QED 
and self-consistent QED were given. Again, these rules are 
valid also in the presence of bound electrons, because the 
bound states do not modify the short-distance or large- 
momentum asymptotics of the Green functions. 

Much work remains to be done to transform QED 
into a practical tool for solid-state or molecular calcula- 
tions. We must investigate more fully the conditions of 
renormalization. We should also be able to calculate the 
energy-momentum tensor from the Green functions, as the 
current was written in terms of the electron propagator. 
Results in this direction were obtained by Engel and Drci- 
zler [|l], but much work remains to be done to reach an 
all-order renormalization of the QED energy. The energy 
is not the only physically interesting property and the in- 
vestigation of the full energy-momentum tensor might be 
interesting and might provide some physically reasonable 
conditions of renormalization. 

In this paper, we have only considered a coherent ex- 
ternal source j^{x). It is also possible to derive Schwinger 
equations for a partially coherent source. The main appli- 
cation of this extension is the effect of temperature on an 
electronic system. 



Finally, it would be interesting for the spectroscopic 
applications to generalize the Schwinger equations to a 
degenerate "vacuum" . When the unperturbed state is de- 
generate, we must consider not only a single matrix ele- 
ment (^\A^(x)\^) but a matrix ((P^A^x^j) . The modi- 
fications induced by the presence of a degenerate vacuum 
will be presented in a forthcoming publication. 
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14 Appendix 

In this appendix, we give the proof of equations that are 
used in the text. 

14.1 Derivative of S%{A) 

We saw in section | that S%(A) = S$(A) + P N with 

0<E n <E F 

For notational convenience, we suppress the argument A 
of the Green functions S. In this equation, the sum is over 
N states. The wavefunctions ip n (x) are solutions of the 
Dirac equation in the presence of A M (x): 

{ihc~f ■ d — mc 2 — e7 • A(x)) ip n (x) = 0. 
Therefore 

(ihcj ■ d — mc 2 — ■ A(x)) P N (x,y) = 0, 

and the adjoint equation, deduced from Eq.(pT[) 

( - iked" - eA»(x))P N (y, x) lfi ~ mc 2 P N (y, x) = 0. 

We perturb the vector potential into ^4p(x) + eV^(x) so 
that Ppf becomes Pjv + e^-P/v, which must be a solution of 
the Dirac equation for the perturbed potential. Keeping 
only the terms linear in e we obtain 

(«fic7 • d — mc 2 — e"f ■ A(x)^SPn(x, y) = ey ■ V(x)Pm(x, y). 

The general solution of this equation is 

6P N = eS%{A) 1 -VP N +R N , 

where Rn is a solution of the unperturbed Dirac equation. 
We repeat the argument for the adjoint equation and we 
obtain 

6P N = e^ °(A) 7 • VP N + eP Nl ■ VS$(A) + R' N , 
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where R' N is a solution of the Dirac and adjoint Dirac 
equations. Finally, we calculate the functional derivative 
by taking V — 8\^8{x — z) so that 

5P N (x,y) A 

+P N (x,z) 1 x S° (z,y) + R / N . 

To determine R' N we require that the perturbation does 
not change the number of bound electrons. In other words 



AN = / dx tr 



7 



)6Pn(x,x) 



0. 



SA x (y) 

To evaluate this, we first calculate 

dxtr[ 7 5 (x, 2 /) 7 A P J v(2/,x)]. 

We use Eq.(|22|) for E F — 0, because of the cyclic property 
of the trace we can bring the last ip m (x) in the front, 
then we use the orthogonality of the solutions of the Dirac 
equation and we find 



dxtr[y , 5 °( 2 ;,y) 7 A P A ,(y,x)] = 6(x« - y u ) 

x tT bPm(yh x ipm(y)]- 

0<E m <E F 



The same equation for the second part (PjvSq) gives the 
same expression, where 9(x° — y°) is replaced by 8(y° — x°). 
Thus 



dxtr[ 7 u (S u (x, y)^P N (y, x) + P N (x, yh A S»(y, x))} 

= Y tr bPm(vh x ^m{y)]- 

0<E m <E F 



Now if we take R' N (x,y; z, A) = Pn(x, z)^ x P^(z,y) we 
obtain 



dxtx[y R' N (x, x; y, A)] = - ^ tJ# m (y) 7 Vm(y)]- 

0<E m <E F 



which compensates exactly for the previous term. Finally, 
since we know that 



5Sg(x,y;A) 



5A\(z 
we have obtained that 



= S °(z,z;.4) 7 A S (z,2,;,4) 



SS N (x, y; A) _ _ ,\ c n 



^(x,z;A) 7 A ^(z,y;A). 



This is a very satisfactory result, which shows that the 
perturbative solution of the Schwinger equation has the 
same form for the vacuum and the charged cases. 



14.2 Integral 

In this section we show the following identity 

F(<p) = F(0) + f&j dx6 ~J^^- ( 51 ) 
This is a functional form of the classical Taylor formula 

f(x) = /(Q) + / d\f'(Xt)t, (52) 
Jo 

where fit) is the derivative of fit). The n-dimensional 
generalization of Eq. (^2|) is 

/(x) = /(0) + jf dA^/(Ax) 



= /(0) + 



\ dAV^/(Ax) 

J 4=1 



(53) 



Eq.(51) can be derived as the above n-dimensional 
case, but wee choose to prove it by an explicit calcula- 
tion which shows that the purpose of the integral over A 
is to change the multiplicity of some terms. We start from 
the Taylor expansion for functional derivatives 



S n FiO) 



5tp(y n ) ■ ■ ■ 5tp(yi) 



x<p(yi) ■ ■■<p(y n )- 



The notation used in this equation means that the func- 
tional derivatives of F are taken at <j> = 0. Therefore, 
using the symmetry of functional derivatives with respect 
to their arguments 



6F(<p) g I 



5ip(x) ^(n-1)! 



dyi • ■ ■ dy n -i 



S n FiO) 



6tp(x)8(p(y n -i) ■ ■'Stp(y 1 ) 



<p(yi) ■ ■■vivn-i)- 



In Eq.(|5l|) the purpose of the integral over A is to re- 
place the factor l/(n — 1)! by the correct factor \/n\. 



SFjXip) _ A"' 1 
8ip{x) ^ 



n=l v ' 



dy\ ■ ■ -ch/ n _i 



(S n F(0) 



-(piyi) ■ ■■<fiy n - 



8cp{x)8(piy n ^ 1 ) ■ ■•5<p{yi) 

The notation SFiXip)/8ipix) has the following meaning. 
The function G(a:, ip) = SFiip) / f 6tp(x) is a function of two 
independent variables (x and ip). Then 8FiX(p)/8(fix) = 
G(x,\<p). 
Therefore, 



8F{Xip) 



o ^ Scp(x) 



If 

<P(%) = —{ I dyi--- dy n -\dx 
8 n F(0) 



n—l 
X 



8ip(x)8cp(y n -i) ■ --Sipiyi) 
XLp( yi ) ■ ■•y>(y n _ 1 )^(x) 
F(<p)-F(Q). 



Christian Brouder: Renormalization of QED in an external field 



21 



Note that an obvious consequence of Eq.(|5l|) is that 
8 



<V0) Jo 
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